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Abstract 

Recently, it has been proposed that the turbulent fluctuations measured in a linear plasma 
device could be described as a superposition of uncorrelated Lorentzian pulses with a narrow 
distribution of durations, which would provide an explanation for the reported quasi-exponential 
power spectra. Here, we study the applicability of this proposal to edge fluctuations in toroidal 
magnetic confinement fusion plasmas. For the purpose of this analysis, we introduce a novel 
wavelet-based pulse detection technique that offers important advantages over existing techniques. 
It allows extracting the properties of individual pulses from the experimental time series, and 
quantifying the distribution of pulse duration and energy, as well as temporal correlations. 

We apply the wavelet technique to edge turbulent fluctuation data from the W7-AS stellarator 
and the JET tokamak, and find that the pulses detected in the data do not have a narrow distri- 
bution of durations and are not uncorrelated. Instead, the distributions are of the power law type, 
exhibiting temporal correlations over scales much longer than the typical pulse duration. These 
results suggest that turbulence in open and closed field line systems may be distinct and cast doubt 
on the proposed ubiquity of exponential power spectra in this context. 
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The understanding of plasma turbulence is of vital importance for the development of 
fusion as an alternative energy source, in view of the fact that it is the dominant mechanism 
behind the unwanted energy loss of confined plasmas, referred to as power degradation. 
Over the last three decades, the complex dynamics of turbulent transport in magnetically 
confined toroidal devices such as tokamaks and stellarators have been probed experimentally 
by examining edge fluctuation data measured with Langmuir probes. It was found that 
edge fluctuations exhibit self-similar properties in certain regimes, leading to power spectra 
characterized by power law sections extending over a meaningful range of frequencies l|-l5|. 

Recently, Maggs et al. jsl have proposed that this type of results should be reinterpreted 
in terms of an uncorrelated superposition of Lorentzian pulses with a narrow distribution 
of pulse durations. The resulting power spectrum would be predominantly exponential, 
which would fit many of the spectra reported in literature [t-^. According to the authors, 
a consequence of this proposal is that the underlying dynamics would be low- dimensional 
quasi-deterministic chaos. This view appears to deny the importance of the mutual in- 
teraction between confinement and turbulence, widely considered a vital ingredient in the 
understanding of power degradation in fusion plasmas, and contrasts with the hypothesis of 
self-organized critical dynamics proposed in the past [lo[ 

In this Letter, we study this proposal using edge fluctuation data from a tokamak and a 
stellarator. As the pulse distribution is central to the issue at hand, we introduce a novel 
wavelet decomposition technique that allows recovering its complete statistical properties. 
The analysis clarifies that there are significant differences between turbulence in linear plas- 
mas and toroidal fusion plasmas, at least in the operational regime examined (L-mode). 

We hold that the power spectrum of a signal is not determined by the spectrum of the 
individual pulse, unless it is composed of uncorrelated pulses with identical duration. Any 
experimental signal s(t) can be approximated by a sum of pulses: 



Here, L{t) is the basic pulse shape, and t,, i = 1, . . . , is a set of times at which pulses 
of amplitude and duration occur. For discretely sampled data, the absolute difference 
between the signal s and the sum s' can be made arbitrarily small, provided L{t) satisfies 
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some minimal requirements. For fixed = r, Eq. ([T]) can be rewritten as 
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such that s'{t) is a convolution of the pulse shape L(t) and a function (in square brackets) 
specifying the amplitude and time of the pulses. Consequently, the power spectrum of s' is 
the product of spectrum of this function and the spectrum of the pulse shape. 

It is useful to choose L(t) as similar as possible to the typical shape of actual observed 
pulses in the data series - which can in principle be determined using techniques such as 



conditional averag: 



ing ^,^1 



14| . Following Maggs et al. {g], one could use L{t) = l/7r(l+t^ 



(Lorentzian). However, an alternative pulse shape such as G{t) = 1/(1 + e'*') provides an 
equally good fit to the experimental pulse shape obtained from our data using conditional 
averaging. While the spectrum of L{t) is exponential {L{u) oc e"'"^'), the spectrum of G 
is power-law (lim^_j.oo (5(0;) oc uj~^). This seems to suggest that the pulse shape does not 
predetermine the shape of the experimental power spectrum. 

Wavelet analysis. Next, we introduce a wavelet technique to obtain the statistics of and 
the correlation between pulses of an arbitrary time signal, with potential application in wide 
variety of fields. We construct a ('mother') wavelet by combining two pulses with opposite 
sign: ipi(t) = L(t) — ^L (|). By analogy with the well-known Difference Of Gaussians (DOG) 
wavelet, we call this new wavelet the Difference Of Lorentzians (DOL) wavelet. Its Fourier 
transform is i>i{co) = (exp(— |a;|) — exp(— 2|a;|))/-\/27r (essentially, exponential decay). For 



this wavelet, the admissibility constant 



15| = I \ipi{uj)\'^/\uj\ duj satisfies < < 



oo, leading to a proper wavelet transform. This makes the full power of continuous wavelet 
analysis available to the issue at hand 16|. We define the wavelet by iprif) = '^i^i (7)) 
such that its L^-norm does not depend on r. The wavelet transform and its inverse are 



a(t,r) = / s{t')^ljr{t - t')dt\ 
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This allows decomposing any (square integrable) signal s{t) into sums of Lorentzians with 



arbitrary precision 



15|] - in particular, any such signal s{t) with any given spectrum. Conse- 



quently, a signal composed of a sum of Lorentzians, Eq. ([T]), can have an arbitrary spectrum. 

The conditions to obtain an exponential power spectrum from the pulse distribution 
are easily gleaned from Eq. (jl]). By setting all wavelet coefficients a to zero except those 
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corresponding to a specific scale r (formally, a{t,r) = a(t,T)S{T — t)), the expression for 
the reconstructed signal s' becomes a convolution, and the power spectrum Pg' = can 
be written: 

PAuj)=Paiuj)-P4^), (5) 

Thus, provided the spectrum of a{t,T) is essentially fiat (i.e., a is random, uncorrelated) , the 
spectrum of s' will have the same shape as the wavelet spectrum (i.e., essentially exponen- 
tial). Note that a similar conclusion can be obtained directly from Eq. ([2]), without invoking 
wavelets. Thus, it is not sufficient to require only that the signal consists of Lorentzian 
pulses with a fixed duration in order to produce an exponential spectrum. This is easily 
demonstrated by generating a signal based on the following simple algorithm: (1) Choose 
a number of pulse amplitudes Oj, i = 1, ...,N. (2) Choose — 1 time intervals At,,. (3) 
Construct the pulse time array from ti = + Atj_i. (4) Compute the signal s{t) from 
Eq. ([1]), using = r. The resulting signal is a sum of fixed-duration Lorentzian pulses. 

By way of example, we set = 5000 and r = 3. For Atj we take a uniformly distributed 
random variable in the range [0, 10), and for Cj a fractional Gaussian noise (fCn) with a 



certain degree of persistence (Hurst parameter H = 0.8) 17|]. The spectrum of this artificial 
signal is shown in Fig. [H It clearly follows a power-law for low frequencies (extending over 
nearly 3 orders of magnitude of /), while at high frequencies it falls of exponentially (as 
expected). While we do not pretend this simulation to be a model of any physical system, 
it certainly shows that a signal consisting of a sum of fixed duration Lorentzian pulses does 
not possess an exponential spectrum unless the pulses are uncorrelated in time. 

Analysis of stellarator/tokamak turbulent fluctuation data. The wavelet trans- 
form, Eq. ([3]), offers an efficient technique for identifying pulses in experimental data s{t). 
For this purpose, the transform a{t,T) is computed and the local maxima of |a(t,r)| are 



found [l8|. The set of values Cmax, ^max, Tmax indicate, respectively, the amplitude, time of 
occurrence, and duration of pulses in the data series. This technique is capable of identify- 
ing all individual pulses, regardless of their scale or amplitude, unlike, e.g., conditional or 
threshold techniques - assuming, of course, that the wavelet shape ipi adequately matches 
the experimental pulse shape and that pulses do not overlap. 

To demonstrate the technique and test the hypothesis of Maggs et al. Fig. |2] shows a 
wavelet spectrum for Langmuir probe data (the ion saturation current, Jsat) obtained in the 
plasma edge of W7-AS, discharge 35425, 430 < t < 455 ms These data are obtained 
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FIG. 1: Spectrum for artificial data according to Eq. ([T]), constructed using the algorithm described 
in the text, showing that even a signal built from fixed-duration Lorentzian pulses is not always 
completely exponential. 

just inside the last closed magnetic surface in an L-mode discharge. 

Fig. |3] shows the probability distributions of the (absolute) amplitude ^(ainax), the 
pulse duration p^r^a.^), and the waiting time between pulses p(Atinax)- To calculate 
the waiting time distribution, only energetic pulses are considered (i.e., pulses satisfying 
E = (amax/Tmax)^ > Eg, whcrc Eg IS an appropriate threshold 19|]). All show a power law 
structure. It is clear that neither the assumption of a narrow pulse duration distribution, 
nor that of an uncorrected pulse train appear justified (the latter would imply a Poisson 
waiting time distribution). 



ts 



or /sat measurements in the edge of a JET plasma (discharge 



20| . We draw attention to the significant length of the observed 



Fig. H] shows similar resul 
54008, 62.52 < t < 62.58 s) 
power laws. 

Discussion. In previous work, the shape of experimental spectra in the edge of fusion 
devices was studied, and power laws were detected over a range of frequencies in certain 
conditions (typically, measurements obtained in L-mode discharges at positions either inside 



or outside of, but very close to, the last closed surface) 



l|-|5|. Many factors may affect the 
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FIG. 2: Top: a short time section (2 ms) of /gat data. Bottom: wavelet transform |a(i,r)| (blue: 
small amplitude; red: large amplitude). The vertical axis shows / = 1/r; only the lowest frequen- 
cies are shown (/ < 100 kHz), although the spectrum is computed up to the Nyquist frequency. 
White spots indicate the detected pulses (local maxima of |a(t,r)|). 

shape of the spectra, such as the Doppler effect due to plasma rotation, the presence of 
modes, and the position of the probe in the plasma, making its study rather non-trivial. 

The analysis presented here, based on a novel technique that is more powerful than 
techniques used in the past, provides deeper insight into the complex nature of the dynam- 
ics in toroidal magnetic confinement fusion plasmas: the distribution of pulse amplitudes 
and durations follows a power law, and long-range temporal correlations (memory effects) 
are significant. Memory effects refer to persistence over times much longer than the tur- 
bulence autocorrelation time, presumably caused by the self-organization of profiles and 



turbulence 



21| . The existence of such memory effects is also supported convincingly by a 



number of studies invo 
cross correlation 



[l|,|22[J25 



ving techniques such as the rescaled range (Hurst parameter 



25| . the structure function |26|, |27|, and quiet time analysis 



28 



and 



29|. 



In these studies, temporal correlations are reported for an ample set of toroidal fusion de- 
vices and it was concluded that memory effects can exceed the local decorrelation time 
by up to three orders of magnitude. 

Conclusions. In this work, we have examined the claim put forward in that pulse 
shapes and spectra of fluctuations in the edge of fusion plasmas might be explained by a 
sequence of Lorentzian pulses, associated with exponential spectra. To address this issue, a 
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FIG. 3: Analysis of the pulse distribution obtained from wavelet analysis (W7-AS). Top: pulse 
amplitude p(amax)- Center: pulse duration p(Tmax)- Bottom: waiting time between energetic pulses 

new pulse detection technique was developed based on a wavelet closely matching the typical 
observed pulse shapes. This technique may find application in other fields where complex 
dynamics are thought to play an important role. 

The analysis has revealed that exponential spectra are only obtained from Lorentzian 
pulses if: (i) the distribution of pulse durations is narrow and (ii) the temporal distribution 
of pulses has a flat spectrum (implying no temporal correlations). The distribution of 
pulse amplitudes and durations calculated for actual fusion plasma e 
temporal correlations are strong, in accordance with earlier work [l| 
and invalidating both these requirements. 
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FIG. 4: Analysis of the pulse distribution obtained from wavelet analysis (JET). Top: pulse 
amplitude p(amax)- Center: pulse duration p(Tmax)- Bottom: waiting time between energetic 

pulses p(Atmax)- 

Thus, we conclude that the dynamics of turbulence in the edge of magnetically confined 
toroidal plasmas are often quite different from and richer than in linear plasmas. In linear 
devices, it seems likely that parallel losses along the field lines to end plates are so fast 
that they limit the range of times over which temporal correlations can be established, 
whereas in toroidal plasmas the existence of closed magnetic surfaces allows phenomena with 
longer time scales to manifest themselves. In this situation, the interaction between profiles 
and turbulence (plasma self-organization) may then be a source of long-range temporal 



correlations 2l|, which may require conditions of good confinement and strong gradients. 



Although the results reported here seem to be typical for tokamaks and stellarators, it 
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cannot be excluded that there are specific regimes and/or regions in toroidal plasmas that 
might exhibit dynamics more similar to the results reported by Maggs et al at LAPD-U. 
For instance, the decorrelation between pulses should be faster in a region of strong poloidal 
sheared fiows or in the far scrape-off layer, where magnetic field lines are connected to the 
first wall. 
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